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Abstract 

A solution of linearized Einstein field equations in vacuum is given and 
discussed. First it is shown that, computing from our particular metric the 
linearized connections, the linearized Riemann tensor and the linearized 
Ricci tensor, the linearized Ricci tensor results equal to zero. Then the 
effect on test masses of our solution, which is a gravitational wave, is 
discussed. In our solution test masses have an apparent motion in the 
direction of propagation of the wave, while in the transverse direction 
they appear at rest. In this way it is possible to think that gravitational 
waves would be longitudinal waves, but, from careful investigation of this 
solution, it is shown that the tidal forces associated with gravitational 
waves act along the directions orthogonal to the direction of propagation 
of waves. The computation is first made in the long wavelengths approx- 
imation (wavelength much larger than the linear distances between test 
masses), then the analysis is generalized to all gravitational waves. 

In the last sections of this paper it is shown that the frequency depen- 
dent angular pattern of interferometers can be obtained from our solution 
and the total signal seen from an interferometer for the stochastic back- 
ground of gravitational waves is computed. 

PACS numbers: 04.80.Nn, 04.30.Nk, 04.30.-w 
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1 Introduction 



The design and construction of a number of sensitive detectors for gravitational 
waves (GWs) is underway today. Some laser interferometers like the Virgo detec- 
tor, being built in Cascina, near Pisa by a joint Italian- French collaboration, the 
GEO 600 detector, being built in Hanover, Germany by a joint Anglo-Germany 
collaboration, the two LIGO detectors, being built in the United States (one in 
Hanford, Washington and the other in Livingston, Louisiana) by a joint Caltech- 
Mit collaboration, and the TAMA 300 detector, being built near Tokyo, Japan, 
are going to become operative in the next years. Many bar detectors are cur- 
rently in operation too, and several interferometers and bars are in a phase 
of planning and proposal stages (for the current status of gravitational waves 
experiments see [H[2])- 

The results of these detectors will have a fundamental impact on astrophysics 
and gravitation physics. There will be lots of experimental data to be analyzed, 
and theorists will be forced to interact with lots of experiments and data analysts 
to extract the physics from the data stream. 

Detectors for GWs will also be important to verify that GWs only change 
distances perpendicular to their direction of propagation and to confirm or ruling 
out the physical consistency of General Relativity or of any other theory of 
gravitation [HUH [5]. 

The response of interferometers to GWs has been analyzed in lots of works 
in literature especially in the so called transverse-traceless (TT) solution of 
linearized Einstein field equations (LEFEs) in vacuum. In this paper a different 
solution of LEFEs in vacuum is given and discussed. First it is shown that, by 
computing from our particular metric the linearized connections, the linearized 
Riemann tensor and the linearized Ricci tensor (for details about linearized 
quantities see [6]), the linearized Ricci tensor results equal to zero (i.e. our 
metric is solution of LEFEs). Then the effect on test masses of this solution, 
which is a GW, is discussed. In our solution test masses have an apparent 
motion in the direction of propagation of the wave, while in the transverse 
direction they appear at rest. In this way it is possible to think that GWs 
would be longitudinal waves, but, from careful investigation of our solution, 
it is shown that the tidal forces associated with GWs act along the directions 
orthogonal to the direction of propagation of waves. The computation is first 
made in the long wavelengths approximation (wavelength much larger than the 
linear distances between test masses i.e. under this assumption the amplitude 
of the GW, h, can be considered "frozen" at a value ho), then the analysis is 
extended to all GWs using a generalization to our solution of the analysis which 
has been used for scalar waves in [5] and for tensorial waves in [7]. 

At the end of this paper it is shown that, from our solution, the angular 
pattern of interferometers can be obtained with a further generalization of the 
analysis of [H[7] and the total signal seen from an interferometer for the stochas- 
tic background of GWs is computed. 
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2 A solution of linearized Einstein field equations 
in vacuum. 



Starting from the metric (we work with c = 1 and h = 1 in this paper) 

ds 2 = [1 + h(t - z)](-dt 2 + dx 2 + dz 2 ) + [1 - h(t - z)]dy 2 , (1) 

where h(t — z) <C 1 is a perturbation of the flat Lorentz-Minkowski background, 
let us call R^ypa and R^ v the linearized quantity which correspond to R^ vpa 
and Rf_ iv . Computing them to first order in a general perturbation h pu , it is 
simple to see that the linearized Riemann tensor is given by [6]: 

Then, because using the metric {1} it is 

hoo = /in = ^33 = h, h 2 2 = -h (3) 

and the other components are equal to zero, it is simple to show that the 
linearized field equations 

= 0, (4) 

are satisfed. 

Thus it is possible to say that the metric HJ is a solution of LEFEs in 
vacuum which describes a gravitational wave propagating in the z+ direction. In 
particular the metric ([I]) is a solution for the + polarization of the gravitational 
wave, in fact equation ([I]) can be obteined directly from the + polarization of 
the TT solution [HI [9] with the substitution 

x — > x 



y 

(5) 

z+\H(t-z) 
t - \H{t - z), 



where 



/t — z 
h{v)dv. (6) 
-oo 

Now we discuss the effect on test masses of this solution. Equation ([I]) can 
be rewritten as 

&-&-& = Th + l T^& <r> 
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where r is the proper time of the test masses. 

To derive the geodesic equation of motion for test masses (i.e. the beam- 
splitter and the mirrors of an interferometer) the equation 



dui 1 dg k i k i 



dr 2 dx 
can be used [8]. 

Thus, from the metric ([J) we obtain 



0, (8) 



d-.r 
Tr 1 

d 2 y 
It 1 

d 2 t 



(9) 



The first and the second of eqs. ((9]) can be immediately integrated obtaining 

Ci = const. (10) 



(7.r 
d7 



^ = C7 2 = const. (11) 
ar 

Assuming that test masses are at rest initially we get C\ = C-2 = 0. Thus, 
even if the GW arrives at test masses, there is no motion of test masses within 
the x — y plane. This could be directly understood from eq. {J) because the 
absence of the x and of the y dependences in the metric implies that test masses 
momentum in these directions (i.e. C\ and Ci respectively) is conserved. This 
results, for example, from the fact that in this case the x and y coordinates 
do not esplicity enter in the Hamilton- Jacobi equation for a test mass in a 
gravitational field (see [8] for details). 

Now eq. (0 reads 

In this way, eqs. ((9J begin 



Pk = 

(13) 

d 2 t 1 d t (l+h) 

It 1 2 (1 + h)' 2 

d 2 z 1 d;(l+h) 

dr 1 ~~ 2 (1+h) 2 ■ 
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Now it will be shown that, in presence of a GW, there will be motion of 
test masses in the z direction which is the direction of the propagating wave. 
An analysis of eqs. fl3|l shows that, to simplify equations, the retarded and 
advanced time coordinates (v, w) can be introduced: 

v = t — z 

(14) 

w = t + z. 



From the third and the fourth of eqs. (|T3[) we have 
d dv d w [l + h(v)] 



dr dr [1 + h(v)f 



0. (15) 



Thus we obtain 



dv , . 

-=a, (16) 



where a is an integration constant. From eqs. (fl~2f and (fl6j) . we get 

dw 8 
dr ~ 1 + h 

where 8 = — , and 



(17) 



r = 8v + 7, (18) 

where the integration constant 7 correspondes simply to the retarded time 
coordinate translation v = t — z. Thus it can be put equal to zero without loss of 
generality. Now let us see what is the meaning of the other integration constant 
(3. From eqs. ((T6]l and (fT7|) the equation for z can be written: 

dz 1 8 2 

i'r/rrk-v- <19 » 



When it is h = (i.e. before the GW arrives at the test masses) eq. (|19l) 
becomes 

Tr " ^ " <20) 
But this is exactly the initial velocity of the test mass, thus we have to choose 
8 = 1 because we suppose that test masses are at rest initially. This also imply 
a = 1. 

To find the motion of a test mass in the z direction from eq. lfT8|l we have 
dr = dv, while from eq. (fl7f we have dw = Because it is z = we 

obtain 

^ = \i^r h -dv), (21) 

which can be integrated as 
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(22) 

1 rt — z h(v) , 

= z o-2 Loo irm dv > 

where zq is the initial position of the test mass. Now the displacement of 
the test mass in the z direction can be written as 

a 1 ft-zo-Az h(v) , 

Az = Z — Zn = — t ; , i,/ \ UV 

u 2 J — oo l+h{v) 

(23) 

1 ft-z t±l_ dv _ 

2 J — oo l+h{v) 

Our results can be also rewritten in function of the time coordinate t: 



x(t) = x 

y(t) = yo 

(24) 

z(t) = 2 ~ \ d ( v ) 

r(t) = t-z(t). 

Now let us reasume what happens in our solution: it has been shown that in 
the x — y plane an inertial test mass initially at rest remains at rest throughout 
the entire passage of the GW, while in the z direction an inertial test mass 
initially at rest has a motion during the passage of the GW. Thus it could 
appear that the solution ((J) for GWs has a longitudinal effect and does not 
have a transversal one, but the situation is different as it will be shown in the 
following analysis. 



3 Analysis in the long wavelenghts approxima- 
tion 

We have to clarify the use of words " at rest" : we want to mean that the 
coordinates of test masses do not change in the presence of the GW in the x — y 
plane [5j [6l [3 [9], but it will be shown that the proper distance between the 
beam-splitter and the mirror of our interferometer changes even though their 
coordinates remain the same. On the other hand, it will be also shown that the 
proper distance between the beam-splitter and the mirror of our interferometer 
does not change in the z direction even if their coordinates change in the solution 
©• 

A good way to analyze variations in the proper distance (time) is by means 
of "bouncing photons" : a photon can be launched from the beam-splitter to be 
bounced back by the mirror (see 0[7] and figure 1). 

In this section we only deal with the case in which the frequency / of the 
GW is much smaller than ^- = where 2Tq = 2Lq is the total round-trip 



C 



Y 



Photon 




Figure 1: photons can be launched from the beam-splitter to be bounced back 
by the mirror 



time of the photon in absence of the GW. The analysis will be generalized to 
all frequencies in the next section. 

We assume that test masses are located along the x axis and the z axis of 
the coordinate system. In this case the y direction can be neglected because the 
absence of the y dependence in the metric (TJ implies that photon momentum 
in this direction is conserved [5j \7\ [8] and the interval can be rewritten in the 
form 

ds 2 = [1 + h(t - z)](-dt 2 + dx 2 + dz 2 ). (25) 

Let us start by considering the interval for a photon which propagates in the 
X axis. We note that photon momentum in the z direction is not conserved, for 
the z dependence in eq. ([1]) (HO [8]. Thus photons launched in the x axis will 
deflect out of this axis. But here this effect can be neglected because the photon 
deflection into the z direction will be at most of order h (see [7]). Then, to first 
order in h, the dz 2 term can be neglected. Thus, eq. lf25|) can be rewritten as 

ds 2 = (l + h)(-dt 2 ) + (l + h)dx 2 . (26) 
The condition for null geodesies (ds 2 = 0) for photons gives 

^ Xp j° t0n = ±1 => Xphoton = const ± t. (27) 

In our solution the x coordinates of the beam-splitter and the mirrors are 
unaffected by the passage of the GW (see the first of eqs. 1(23)1 ), then, from eq. 
|27|) . it is possible to see that the interval, in coordinate time t, that the photon 
takes for run one round trip in the x arm of the interferometer is 

T = 2L (28) 
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(i.e. the photon leaves the beam-splitter at t = and returns a t = T). But 
this quantity is not invariant under coordinate transformations, and we have to 
work in terms of the beam-splitter proper time which misures the physical lenght 
of the arms. In this way, calling r(t) and Zb{t) the proper time and z coordinate 
of the beam-splitter at time coordinate t with initial condition z&(— oo) = 0, 
from eqs. (|24|) it is possible to obtain 

z b\P) — 2J-00 i+h(v) au 

(29) 

r(t) = t +i f^ 26(t) T^rdw. 

Thus, calling t x the proper time interval that the photon takes to run a 
round-trip in the x arm, we have 

r x = r(T) - r(0) = T + \ ^dv ~ 

~T+i/i [T + z b (0)-z b (T)]~ (30) 

~2i (l + |M- 

In the above computation eq. ([28)) has been used and, considering only the 
first order in h with h <C 1, the field ft. has been also considered "frozen" at a 
fixed value ho. We note that z&(0) — Zb(T) is second order in ho- 

Now let us consider the z direction: the x direction can be neglected be- 
cause the absence of the x dependence in the metric HJ implies that photon 
momentum in this direction is conserved [21E1H]. From eq. (|25|) it is now: 

ds 2 = (l + h)(-dt 2 ) + {l + h)dz 2 , (31) 
and the condition for null geodesies (ds 2 = 0) for photons gives 

^~ P ^° t0n = ±1 => z photon = const ± t. (32) 

Supposing that the photon leaves the beam splitter at t = let us ask: how 
much time does the photon need to arrive at the mirror in the z axis? Calling 
T\ this time we need the condition 

z b (0)+T 1 = z m (T 1 ), (33) 

where z m (t) is the z coordinate of the mirror in the z axis at coordinate time 
t with z m (— 00) = Lq. In the same way, when returning from the mirror the 
photon arrives again at the beam-splitter at t = T z = T\ + T 2 , then 

z m {T 1 )-T 2 = z b {T z ). (34) 
Subtracting eq. (|34| from eq. ([33f we obtain 

T z =Ti+T 2 = [z m {T x ) - z b (0)} + MTi) - Zb(T z )]. (35) 
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It is known from eq. (|24|) that the equations of motion for z& and z m are: 



*w = -i/-r (i) TTgy^ 

and, substituing them in eq. lj35|). we get 



(36) 



2i_ 2b(0) l + h(v) 2j Tz _ Zb[T ^ l + h(v) y > 

From eq. (f33|l we see that the first integral in eq. |37|) is zero. The second 
integral is simple to compute if the G W is considered frozen at a value ho . To 
first order in this value it is 



~\ StI-zT{£) T+$v) dv - -l h o[Ti - z m {Ti) -T z + z b (T z )} ~ 

(38) 

— — 5^o(-^o — L — 2L ) = +i/io2L - 
In this way eq. (|37|l becomes 

T, = (1 + ^o)2i - (39) 

Then, calling r z the proper time interval that the photon takes to run a 
round-trip in the z arm, with the same way of thinking which leaded to eq. 
(|30]). we obtain 

r z = r{T z ) r(0) = T z + \ ^dv * 

~T z + ±h o {T z + z b (0)-z b (T)}~ (40) 
~T z {l-\h Q )~2L Q . 



Thus, from eqs. (|30|l and (|40j) it is shown that there is a variation of the 
proper distance in the x direction (transversal effect of the G W) , while there is 
not a variation of the proper distance in the z direction (no longitudinal effect). 

4 Generalization of the analysis 

Now the previous result will be generalized to all the frequencies, with an anal- 
ysis that, with a trasform of the time coordinate to the proper time, generalizes 
to our solution the analysis of [7] , where the analysis was made using the TT 
solution of LEFEs in vacuum. In this way it will also be obtained the response 
function of the interferometer for our solution. 
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Let us start with the x arm of the interferometer. The condition of null 
geodesic (|27|) can also be rewrite in this way: 

dt 2 = dx 2 . (41) 

Thus, in this case, the analysis of [7] cannot be used starting directly from 
the condition of null geodesic. In fact the metric f26|) is different from the 
metric of eq. (3) in [7]. In [7] the author used the condition of null geodesic to 
obtain the coordinate velocity of the photon which was used for calculations of 
the photon propagation times between the test masses (eq. (4) in [7]). But in 
equation l|4"l"j) appears that the coordinate velocity of the photon is equal to the 
speed of light in our solution. Then let us ask which is the important difference 
between our metric (Q]) and the TT metric analyzed in 0. The answer is that 
the TT metric defines a "synchrony coordinate system", a coordinate system in 
which the time coordinate t is exactly the proper time (about the synchrony 
coordinate system see Cap. (9) of [8] and Section 4 of [5]). In the solution (TJ 
t is only a time coordinate. We know that the rate dr of the proper time is 
related to the rate dt of the time coordinate from O [8] 

dr 2 =g m dt 2 . (42) 

Only with the aid of the time transform lj42j) the analysis of [7] can be applied 
to our solution. 

From eq. ((26j) it is goo = (1 + h). Then, using eq. (JUJ, we obtain 

dr 2 = {I + h)dx 2 , (43) 

which gives 

dr = ±(1 + h)?dx. (44) 

Now it will be shown that the analysis of [7] works in our case too. This is 
an analysis parallel to the one used in Section 4 of [5]. 

In eqs. ((24)1 it is shown that the coordinates of the beam-splitter Xb = I and 
of the mirror x m = I + Lq do not change under the influence of the GW in our 
solution, thus the proper duration of the forward trip can be found as 

[L +l 

n(t) = J [l + h(t)]vdx. (45) 

To first order in h this integral can be approximated with 

Tl {t)=T + ^J h(t')dx (46) 

where 

f = t- (l + L - x). 
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In the last equation t' is the retardation time (i.e. t is the time at which the 
photon arrives in the position I + Lq, so I + Lq — x = t — t' [5j E]). 
In the same way we have for the proper duration of the return trip, 



T 2 (t)=T + lf h(t')(-dx), (47) 

Z Jl+Lo 



where now 

t' = t-(x-l) 

is the retardation time and 

Tq = Lq 

is the transit proper time of the photon in the absence of the GW, which 
also corresponds to the transit coordinate time of the photon in the presence of 
the GW (see eq. lj4Tj) ). 

Thus the round-trip proper time will be the sum of r 2 (t) and Ti{t — T ). 
Then, to first order in h, the proper duration of the round-trip will be 

TV.t.(t) =n(t-T Q )+T 2 (t). (48) 

Prom eqs. |46| and (|47|) it is immediately shown that deviations of this 
round-trip proper time (i.e. proper distance) from its imperurbated value are 
given by 



1 



L +l 



$ T (t) = o / ~ 2T ° + x-l) + h{t-x + l)]dx. (49) 

2 J i 

Eq. lj49j) generalizes eq. l|30|) which was derived in the low frequencies 
approximation. The signal seen from the arm in the x axis can be also defined 
like 



To 



/ [h{t -2T +x-l) + h{t-x + l)]dx. (50) 
21 o Jl 

Now the analysis will be transled in the frequency domain using the Fourier 
transform of the field h defined by 



h{uj) = / dt h{t)exp{iuJt). (51) 

./ — OO 

With the definition (|5T)l . from eq. (f50|) we have 

^ = T(uj)h(uj), (52) 
where T(w) is the response of the x arm of our interferometer to GWs: 
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z 




Figure 2: the beam splitter and the mirror are located in the direction of the 
incoming GW 



which is computated in lots of works in literature. 

Now let us see what happens in the z coordinate (see figure 2). 

Eq. (|3Tjl and the condition ds 2 = for null geodesies also give 

dz = ±dt. (54) 
But, from the last of eqs. i[24)l the proper time can be written like 



dr(t) =dt- dz, (55) 
and, combining eq. i|54[l with eq. (|55| . we get 

dr(t) =dt^ dt. (56) 

Thus we obtain now 



n(t)=o (57) 



for the forward trip 
and 



12 




2dt = 2T ( 



(58) 



for the return trip. Then 



r(t) = n(t) + T 2 (t) = 2T . 



(59) 



Thus there is not longitudinal effect (i.e. St = SLq = 0). This is a direct 
conseguence of the fact that a GW propagates at the speed of light. In this 
way in the forward trip the photon travels at the same speed of the GW and 
its proper time is equal to zero (eq. (|57|l ). while in the return trip the photon 
travels against the GW and its proper time redoubles (eq. ([58]) ). 

5 The low-frequencies angular pattern of inter- 
ferometers in the TT solution of linearized Ein- 
stein field equations in vacuum 

Very important in a cosmological scenario is the concept of stochastic back- 
ground of gravitational waves (9j [lOl 2U 221 [13] . 

Roughly speaking these are GWs that are produced by a very large number 
of weak, independent and unresolved sources [9l [lOl US [13]. A stochastic back- 
ground of GWs could be the result of processes that took place in a time very 
close to the Planck era [9j [TTJ, [T2], [13] . But it could be due to more recent pro- 
cesses too. An example can be radiation from many unresolved binary systems 
like neutron stars, white dwarfs and black holes. These more recent contribu- 
tions could overwhelm the primordial parts of the background. In this context 
we know that, in any case, the properties of the radiation will be strongly de- 
pendent upon the source. For example, a stochastic background of relic GWs is 
expected to be isotropic, while gravitational radiation derived by binary sistems 
in our galaxy would be highly anisotropic. Thus we have to see the output of our 
detectors before taking a decision between these two possibilities [TUl 221 23] . 

In which sense sources of the stochastic background are "unresolved" can 
be understood making an analogy with optical sources. In the study of an 
optical source, somewhere in the sky, using a telescope with a certain angular 
resolution, details of the source can be resolved if the angular resolution of our 
telescope is smaller than the angular size of the features of the source. In the 
case of Virgo (and similiar experiments like the two LIGO), the angular size 
of the detector pattern is of order 90° . Thus almost any source could make a 
significant contribution to the detector strain for almost any orientation of both 
detector and source, and the sources are unresolved. When lots of sources that 
give a contribution are present, even if they are pointlike, the resulting signal is 
stochastic. 

Thus, in this scenario, it is simple to understand the fundamental importance 
of the angular pattern of a detector. 
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Let us see what happens in a detector. The total output of the antenna S(t) 
is in general of the form 

S(t)=s(t)+n(t), (60) 

where n(t) is the noise and s(t) is the contribution to the output due to the 
gravitational waves. 

For an interferometer with equal arms of lenght L (3 kilometers in the case 
of Virgo) in the u — v plane, it is 

SL u (t) - SL v (t) 
s{t) = - , (61) 

where SL U>V are the displacements produced by gravitational waves. In the 
TT solution of LEFEs in vacuum (this solution is historically called transverse- 
traceless (TT) , because the gravitational waves have a transverse effect and are 
traceless [6]) the total perturbation of a gravitational wave propagating in the 
positive 2* = zf2 direction and with a wave front parallel to the x — y plane [6] 
is given by 

h a(3 (t -z) = h + (t - z)e+ (fl) + h x (t- «)e*,(fi), (62) 
where e+,(Q) and e*p(Q) are the two polarizations [9j HO ) fT H \T % fT3] : 

+ _ At* a/3 _ „i<V,/3 



y a f (63) 



e^ = ^ + iY. (64) 

It is also known that the relation between the output s(t) and the total signal 
of gravitational waves h a b{t) in the TT solution of LEFEs in vacuum has the 

form o nn Ha HE] 

s(t) = D ab h ab (t), (65) 

where D ab is called detector tensor. For an interferometer with arms along 
the wew directions (not necessarly orthogonal) we have 

D ab = ^(v a v b -u a u b ). (66) 

Thus, in the case of the stochastic bacground of gravitational waves, the 
equation 



habit, X > ) = 

= h Ea /_ ~ duj Is* d &~ h A(uJ, Cl) exp iu(t - Cl ■ x > )ef b (fJ) 
can be used for the total signal \M HH [H [13l [H] , where 



(67) 



= cos (f> sin Ox + sin cj> sin By + cos 6z, (68) 
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with 



dQ = d cos 6d(j> (69) 

and 

c^(n)^[e+(n),e a x 6 (n)], (70) 



are the two polarizations l|63p and (|64|) . 

Putting = in our expansion (i.e. the coordinates of the detector are in 
the origin of our system) , we obtain 



s(t) = — V/ duj / dn~h A (uj,n)exp(iujt)D ab e£ b (n). (71) 
2tt ^ J.^ J S 2 

The correspondent equation in the frequency domain is also given by: 

S(uj) = D ab h ab {Lo), (72) 

which can be rewritten as 

s(w)= [ dn~h A (uj,n)D ab e£ b {n). (73) 

Js 2 

The quantity 

F A (n) ee D ab e£ b (Cl) = Tr{De}, (74) 
is called detector pattern. 

Thus, using eqs. ([63]), ([64)) and (J66j) , combined with eq. |74|) it is simple to 
obtain 

F+(n) = -(1 + cos 2 6») cos 20 (75) 

F x (fi) = - cos 6> sin 20 (76) 

These detector patterns for different polarizations of GWs have been ana- 
lyzed in lots of works in literature (see for example [L2 \ [T3" l fT4 l I15|). 

But there is a problem: eqs. <[75j) and (|76| are not the general form of the 
detector patterns, but they are only a good approximation for long wavelengths 
(i.e. the wavelength of the wave is much larger than the linear dimension of 
the interferometer) [15]. In this approximation the detector can be considered 
pointlike when a gravitational wave is arriving (i.e. we can see the wave "frozen" 
at a value ho). 

In the next Section, with the auxilium of our solution of LEFEs in vacuum, 
the exact frequency - dependent expressions for equations (|75j) and ([76]) will be 
derived in the case of an interferometer with perpendicular arms. 
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6 The detector pattern in the general case in our 
solution. 



In the context of the potential detection of stochastic backgrounds of GWs with 
interferometers in this Section it will be generalized to all the wavelengths the 
concept of angular pattern for interferometers, which is well known for the TT 
solution of LEFEs in vacuum and in the assumption that the wavelength of 
the GWs is much larger than the distance between the test masses (i.e. the 
beam-splitter and the mirrors of the interferometer, see Section 5 and refs. 

pa ma nam!). 

We emphasize that, in the analysis of the angular pattern of interferometers, 
a further generalization of the analysis of [5j [7] will be made. In [7] only the 
simplest geometry case was considered in the purely + polarization of the wave 
while in [5] the purely scalar case was analyzed. Here it will be shown that the 
analysis can be generalized for the more general geometry of the + polarization 
and for the x polarization of GWs too. 

It will be computed the variaton of the proper distance that a photon covers 
to make a round-trip from the beam-splitter to the mirror of an interferometer 
[HE] with a coordinates choice that, for the + polarization gives the line element 
(H), while for the x polarization gives the line element 

ds 2 = {-dt 2 + dz 2 )[l + h x (t - z)] + dx 2 + dy 2 + 2h x (t - z)dxdy, (77) 

which can be obteined applying the substitution |(5]) to the TT solution for 
the x polarization. 

It is simple to see that also the metric |77|) satisfes LEFEs in vacuum. 

Now, with a treatment which generalize to the angular dependence the anal- 
ysis of [7] and wich is parallel to the analysis in Section 7 of [5] , the computation 
will be transled in the frequency domain and the general frequency dependent 
angular patterns of interferometers will be derived. 

We start from the + polarization. In this case the interval ((XJ) has to be 
considered. 

But we recall that the arms of our interferometer are in the it and it 
directions, thus, to compute the line element in the it and it directions, a 
spatial rotation of our coordinates has to be made [5] : 

u = — x cos 9 cos 4> + y sin0 + z sin 6* cos <j> 

v = — x cos 9 sin <j> — y cos <j> + z sin 6 sin cf> (78) 

w = x sin 9 -\- z cos 9 , 

or, in terms of the x, y, z frame: 
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Figure 3: a GW incoming from an arbitrary direction 



x = —u cos 9 cos 4> — v cos 9 sin (f> + w sin 9 

y = u sin — v cos <f> (79) 

z = usin#cos0 + v sin # sin </> + wcos9. 

In this way the GW is propagating from an arbitrary direction ~r* to the 
interferometer (see figure 3). 

The metric tensor transforms like [8] : 

n ik _ dx% dx Jim, (m \ 



Using eq. 1(78)1 . eq. |79)) and eq. (|8Q|) . in the new rotated frame, the line 
element {1} in the T? direction becomes (here we can neglect the v and w 
directions because we will use bouncing photons and the photon deflection into 
the v and w directions will be at most of order h, then, to first order in h, we 
can neglect the dv 2 and dw 2 terms): 



ds 2 = [1 + (cos 2 6» cos 2 </> - sin 2 (f>)h + (t - u sin 9 cos (j))](du 2 - dt 2 ). (81) 

We emphasize that, in the line element f81)) a spatial dependence and an 
angular dependence are present in the u direction, differently from the line 
element of eq. 27 of [7]) where, because of the simplest geometry there is a 
purely time dependence. 

The condition for null geodesies (ds 2 = 0) in eq. f8l)) gives the coordinate 
velocity of the photon: 

du 2 =dt 2 . (82) 
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Thus, also in this case, the analysis of [7] cannot start directly from the 
condition of null geodesic. In fact also the metric lj8lj) is not a "synchrony 
coordinate system". Thus also in the coordinates (|8"T|) t is only a time coordinate 
(i.e. the rate dr of the proper time is related to the rate dt of the time coordinate 
from eq. (|4"2)l ). Then in this case a generalization of our analysis of Section 4 
can be applied like in Section 7 of [5]. 

From eq. ([8T]) it is 



5oo = [1 + (cos 9 cos 4> - sin 2 (j>)h + (t - -usinflcos^)]. (83) 



Then, by using eq. (|82|) we obtain 

dr 2 = [1 + (cos 2 9 cos 2 <f) - sin 2 </>) h + (t-u sin 9 cos (f>)] )du 2 , (84) 
which gives 

dr = ±[1 + (cos 2 9 cos 2 - sin 2 (f>) h + (t - u sin cos 0)] i du. (85) 

Now it will be shown that the analysis of [7] works in this case too. 

We put the beam splitter in the origin of the new coordinate system (i.e. 
Ub = 0, Vb = 0, Wb = 0) [5]. From eqs. (f24j) it is known that an inertial test mass 
initially at rest in the x — y plane in our coordinates, remains at rest throughout 
the entire passage of the GW. Eqs. |2~4|) also show that the coordinates of the 
beam-splitter and of the mirror change under the influence of the GW in the z 
direction, but this fact does not influence the total variation of the round trip 
proper time of the photon (eq. lf59j) ). Then, in the computation of the variation 
of the proper distance in our coordinates, the coordinates of the beam-splitter 
Ub = and of the mirror u m = L can be considered fixed even in the u — v 
plane, because the rotation (|78| does not change the situation. Thus the proper 
duration of the forward trip can be found as 



n(t) 

with 



/ [1 + (cos 2 9cos 2 cj) - sin 2 cj))h + (t - u sin 9 cos cj))}^du. (86) 
Jo 



t' = t— (L — u). 



In the last equation t' is the retardation time (see Section 4 and refs. [11 [7]). 
To first order in h + this integral can be approximated with 

, , „ cos 2 9 cos 2 6 — sin 2 6 f L , , , , . . ,. , . . 

r 1 (t)=T+ ^- -/ h+(t' -u sin 9 cos <j))du, (87) 

where 

T = L 



o 
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is the transit time of the photon in the absence of the GW. Similiary, the 
duration of the return trip will be 



. , ^ cos 2 cos 2 6 — sin 2 6 f , , , , . „ ,., , . 
r 2 (t)=TH -/ /i + (i'-usin(9cos0)(-du), 

though now the retardation time is 

t' = t- (u- I). 

The round-trip time will be the sum of T2(t) and n [i — 72(t)], where T2(t) is 
the coordinate time which corresponds to T2(i). The latter can be approximated 
by T i(t—T) because the difference between the exact and the approximate values 
is second order in h + . Thus, to first order in h + , the duration of the round-trip 
will be 

T r . t .{t)=T 1 {t-T) + T 2 (t). (89) 

Using eqs. lj87j) and (|88|) it appears immediatly that deviations of this round- 
trip time (i.e. proper distance) from its imperurbated value are given by 



ST (Q = cos 2 9 cos 2 0-sin 2 r£ + (f _ 2T + - Sill 6 COS 0)) + 

(90) 

+/i + (t — u(l + sin0cos(/>))]du. 

Now, with the auxilium of the Fourier transform of the + polarization of the 
field, defined by 

/oo 
dth + (t) exp(iujt) (91) 
-co 

we obtain, in the frequency domain: 

Sf(t) = (cos 2 6 cos 2 4> - sin 2 <jj)H. + (w, 0, (a;) (92) 

where 

H+(w,6,(t>) = -i+^p(2^l) 

" V 1 ' 22lj(1 — sin o cos (p) 

(93) 

i sin 9 cos 0(( l+cxp(2icjL) — 2 cxp icjL(l+sin cos 0))) 
' 2iw(l-sin :i Scos 2 0) ' 

and we immediately see that H+(uj, 9,4>) — > L when a; — ► 0. 
Thus, the total response function of the arm of the interferometer in the 
direction to the + component of the GW is: 

% (cos 2 9 cos 2 6 — sin 2 6) ~ , , „ ,. f , , , . , 
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where 2L = 2T is the round trip time in absence of gravitational waves (note 
that in [7] the Laplace transforms have been used. Here we use the Fourier ones 
because the frequency response functions of the Virgo interferometer for the two 
polarizations of the GW will be designed, like in [5] for the scalar case). 

In the same way the line element (I]) in the if direction becomes: 



ds 2 = [1 + (cos 2 9 sin 2 (/>- cos 2 (f>)h + {t~v sin 9 sin (f))}(dv 2 - dt 2 ), (95) 

and the response function of the v arm of the interferometer to the + polar- 
ization of the GW will be: 

, , . (cos 2 9 sin 2 6 — cos 2 d>) - , . . - , . , , , 

T ^)= 2L -K(uj,0^)h+(Lu) (96) 

where now it is 

— l+cxp(2icjL 



u v ' ' 2ilj(1— sm z 9 sin^ <p) 
sin sin 0((l+cxp(22Ct;L) — 2 cxp iuL(l-\-s'm sin </>))) 



(97) 



with iJ+(w, 9,4>) ^ L when u> — > 0. 

Thus the total response function (i.e. the angular frequency dependent de- 
tector pattern) of an interferometer to the + polarization of the GW is: 

fr+t \ (cos 2 9 cos 2 (h— sin 2 <h) fr / /i j\ , 

H + (u) = i ^-H u (u,6,(p)+ 

(98) 

(cos 2 9 sin 2 <b— cos2 4>) tt ( n ±\ 

2T -H v (uj, 9, (/)) 

that in the low frequencies limit (u — > 0) is in perfect agreement with the 
detector pattern of eq. (|75l) : 



-> 0) = -(1 + cos 2 6») cos 20. (99) 

Now the same analysis can be applied to the x polarization. In this case 
the line element (|77|l has to be considered, and, using eq. l(78|). eq. I|79p and eq. 
ll80l). in the new rotated frame the line element (l77l) in the it direction becomes: 



ds 2 = [1 - 2cos6>cos0sm(/>/i x (£ - u sin 9 cos <p)](du 2 - dt 2 ). (100) 

In this way the response function of the u arm of the interferometer to the 
x polarization of the GW is: 

„„ „ , . — cos 9 cos <b sin d> ~ v , „ , . . 
T* (w) = j-^ ^H* (cj, 9, d>). (101) 

In the analogous way the line element (|77f in the if direction becomes: 
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ds 2 = [1 + 2 cos 6 cos cf) sin (f>h x (t - u sin 9 sin 4>))(dv 2 - dt 2 ) (102) 

and the response function of the v arm of the interferometer to the x polar- 
ization of the GW is: 

™v/ x cos cos oi sin r6 ~ . 

T„ x (w) ^ £fT* (w, 0, 0) (103) 

Thus the detector pattern of an interferometer to the x polarization of the 
GW is: 

ttv , \ — cos 6 cos 6 sine!). ~„. „ ,. ~„. „ ... .„ „„. 

ff"(«) = ^ \ H u («> 0, 0) + * (w, A, </>)] (104) 

that in the low frequencies limit (w — > 0) is in perfect agreement with the 
detector pattern of eq. (|76|) : 

ff x (w -► 0) = - cos (9 sin 2i9. (105) 

Then it has been shown that, with the auxilium of our solution of LEFEs in 
vacuum, a generalization of the analysis in [7] works in the computation of the 
two frequency - dependent detector patterns of interferometers exactly like for 
the scalar waves in [5]. 

In figs. 4 and 5 the absolute value of the total response function of the 
Virgo interferometer (L = 3Km) for the + and x polarizations of gravitational 
waves propagating from the direction 6 = \ and ^ = § are respectively shown. 
From the figures it appears that at high frequencies the absolute value of the 
response function decreases respect to the constant value of the low frequencies 
approximation. 



7 The total output due to the stochastic back- 
ground of gravitational waves 

With the auxilium of the frequency dependent pattern functions the approxi- 
mate eqs. (|73)l and (|7T|) can be generalized with, respectively, 

S(u)= f d(l[H x (Lu)h x +H+(uj)h + } (106) 
Js 2 

and 



1 



+ 00 



s(t) = Y I duj dtt[H x {uj)h x + H+(uj)h + }cxp(iujt). (107) 



21 



H(f) 




Figure 4: the absolute value of the total response function of the Virgo interfer- 
ometer to the + polarization of the gravitational waves for 6 = \ and </>=■§■ 



22 



0.38 


; |H(f) | 






. 375 








. 37 








. 365 








. 36 








. 355 










2000 4000 


6000 8000 


10000 


Figure 5: the absolute value of the total response 
ometer to the x polarization of the gravitational 


function of the Virgo interfer- 
waves for 9 = j and <fi = f • 





f 



23 



8 Conclusions 



A solution of linearized Einstein field equations in vacuum has been given and 
discussed. First it has been shown that, if from the metric (U) the linearized 
connections, the linearized Riemann tensor and the linearized Ricci tensor are 
computed, the linearized Ricci tensor is equal to zero (the metric |TJ) is a so- 
lution of LEFEs in vacuum) . Then the effect on test masses of the particular 
solution {J}, which is a gravitational wave has been discussed. It has been 
shown that in our solution test masses have an apparent motion in the direction 
of propagation of the wave, while in the transverse direction they appear at 
rest. In this way it could appear that GWs would be longitudinal waves, but, 
from careful investigation of the solution ([1]) , it has been found that the tidal 
forces associated with GWs act along the directions orthogonal to the direc- 
tion of propagation of waves. The computation has first been made in the long 
wavelenghts approximation (wavelength much larger than the linear dimensions 
of the interferometer) , then the analysis has been applied to all GWs using a 
generalization to our solution of the analysis of [U [7] . 

After this, in the context of the potential detection of stochastic backgrounds 
of gravitational waves with interferometers, which is very important in a cos- 
mological scenario because a stochastic background of GWs could be the result 
of processes that took place in a time very close to the Planck era, our solu- 
tion for LEFEs in vacuum has been used to generalize to all the frequencies 
the concept of detector pattern of interferometers, which was well known in the 
assumption that the wavelength of the GWs is much larger than the distance 
between the test masses (i.e. the beam-splitter and the mirrors of the inter- 
ferometer). In the low frequencies approximation our results agree with the 
standard detector patterns computed in in lots of works in literature (see for 
example [9l [T2l [131 [141 [T5j) , while, at high frequencies, it has been obteined that 
the absolute value of the total response function of interferometers decreases 
with respect to the constant value of the long wavelenghts approximation. 

We emphasize that, also in the analysis of the angular pattern of interfer- 
ometers, the analysis of 0[7] has been generalized: in [7] only the simplest ge- 
ometry case in the TT solution was considered, while in [5] the analysis worked 
for purely scalar waves. It has been shown that the analysis can be generalized 
to the more general geometry and to our solution for LEFEs in vacuum too. 

At the end of this paper we have computed the total signal seen from an 
interferometer for the stochastic background of gravitational waves. 
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